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Abstract 

We study the polarized Bjorken sum rule at low momentum transfers in the range 0.22 < Q < 1.73 GeV with the 
four-loop N 3 LO expression for the coefficient function CBj(a s ) in the framework of the common QCD perturbation theory 
(PT) and the singularity-free analytic perturbation theory (APT). The analysis of the PT series for CBj(a s ) gives a hint 
[to its asymptotic nature manifesting itself in the region Q < 1 GeV. It relates to the observation that the accuracy of 
both the three- and four-loop PT predictions happens to be at the same 10% level. On the other hand, the usage of the 
[two-loop APT allows one to describe the precise low energy JLab data down to Q ~ 300 MeV and gives a possibility for 
reliable extraction of the higher twist (HT) corrections. At the same time, above Q ~ 700 MeV the APT two-loop order 
with HT is equivalent to the four-loop PT with HT compatible to zero and is adequate to current accuracy of the data. 
PACS: ll.10.Hi, 11.55.Hx, ll.55.Fv, 12.38.Bx, 12.38.Cy 
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1. Introduction 

The higher order perturbative QCD (pQCD) and higher 
twist corrections become very important, in particular, in 
observables of the Deep-Inelastic Scattering (DIS) at low 
momentum transfers Q < 1 GeV. The most precise low- 
energy data from the Jefferson Lab [l|, Q on one of the 
main sources of information about the nucleon structure, 
the Bjorken sum rule (BSR) [3j, are the real challenge to 
the accuracy of the pQCD expansions. In our previous pa- 
pers 0, [H, we studied this issue at the three- loop level. In 
the current paper, we continue this line of investigations 
and explore the effect of the recent calculation Q of the 
four- loop (in a s ) contribution to the BSR. 

The BSR claims that the difference of the proton and 
neutron structure functions integrated over all possible val- 
ues 



[gf(x,Q 2 )-g?(x,Q 2 )]dx, (1) 



of the Bjorken variable x in the limit of large four- momentum 
squared of the exchanged virtual photon, Q 2 — > oo, is 
equal to 5,4/6, with gA = 1.267 ± 0.004 7], the nucleon 
axial charge defined from the neutron /3-decay data. 

The r.h.s. of Eq. ([T]) is given by a sum of two series 
in powers of 1/Q 2 (OPE higher twists corrections) and in 
powers of the QCD running coupling a s (Q 2 ) (pQCD radia- 
tive corrections). Until very recently, the pQCD contribu- 
tion to BSR was known @ up to the third order ~ a 3 . So 
far, the corresponding expression has been used in many 



studies, in particular, for extraction of the a s values at low 
momentum scales [9(. 

One of the actual theoretical subjects is the interplay 
between the higher twists (HT) and higher order pQCD 
corrections at low Q, which has recently been studied in 
Refs. 0] at the three-loop level. There, it was shown that 
the satisfactory description of the data down to Q m in ~ 
Aqcd — 350 MeV can be achieved within the Analytic 
Perturbation Theory (APT), the ghost-free modification 
of pQCD. In the current work we repeat this analysis at 
the four-loop N 3 LO level. 

The APT approach is based on the causality principle 
implemented as the analyticity imperative in the complex 
Q 2 -plane for the QCD coupling a s (Q 2 ) in the form of the 
Kallen-Lehmann spectral representation [h| and on the 
demand of compatibility with linear integral transforma- 
tions [ill ('for an overview on the APT concept and results, 
see Ref. [12|]). It is well-known that in the APT framework, 
the theoretical ambiguity associated with pQCD higher- 
loop corrections is diminished (see Ref. [l3|), and results 
are practically renormalization scheme independent [l4j . 

The four-loop expression for the pQCD contribution to 
the Bjorken sum rule became recently available in Ref. @. 
It gives us a reasonable motivation for a new extended 
QCD analysis of the combined JLab data on T\~ n {Q 2 ) at 
low 0.05 < Q 2 < 3.0 GeV 2 accounting for up to o^-order 
in both the (standard) PT and APT approaches. 

The paper is organized as follows. In Section 2, we 
study the higher loop stability of both the PT and APT 
series and the renormalisation scale dependence of the 
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higher-order PT expansion for the Bjorken sum rule. Sec- 
tion 3 contains the QCD results on extraction of the higher 
twist terms from the experimental data at the four-loop 
level. Summarizing comments are given in the last sec- 
tion. 

2. The perturbative QCD contribution 

Commonly, one represents the Bjorken integral ([T]) as a 
sum of the perturbative and the higher twist contributions 
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The perturbative QCD correction^ Aej(Q 2 ) has a form of 
the power series in the QCD running coupling a s (Q 2 ). At 
the up-to-date four-loop (N 3 LO) level in the massless case 
it looks like 



a b j t (Q 2 ) = E ^(Q 2 ). 



(3) 



fc<4 



Here, the numerical expansion coefficients in the modi- 
fied minimal subtraction (MS) scheme, for three active fla- 
vors, n f = 3, read c x = 1/tt = 0.31831, c 2 = 0.36307 0, 
c 3 = 0.65197 and c 4 = 1.8042 @. Besides, the four- 
loop running coupling a s {Q 2 ) is defined as a solution of 
the Renormalization Group (RG) equation 



dL 



= P(a s ); p(a s ) = £ p k , 



k+2 



(4) 



0<fc<3 



where L = ln(/i 2 /A 2 ) and Pk are the coefficients of the 
/3-function. For our purposes, it is convenient to represent 
the /3-function in the form 

f3(a s ) = -A) &l (1 + ha s + b 2 a 2 + b 3 a 3 s +...), (5) 

with bi — Pi/ Po, the ratios of the /3-function coefficients. 
For three flavors the coefficients are fio = 9/4ir = 0.7162, 

b x = 0.5659, bfs = 0.4530 [l6| and bf^ = 0.6770 0. In 
the current analysis we use the exact solutions of the RG 
equation (g]) in the MS-scheme at the scale fx = Q. 

2.1. Analytic Perturbation Theory 

The moments of the structure functions are analytic 
functions in the complex Q 2 -plane with a cut along the 
negative part of the real axis (see, e.g., Ref. [l8j])- The per- 
turbative representation ([3]) violates these analytic proper- 
ties due to the unphysical singularities of a s (Q 2 ) for Q 2 > 
0. To resolve the issue, we apply the APT method [l(J 
[l^ |. which allows one to combine the RG invariance with 
proper analytical properties of the RG-invariant coupling 



lr This correction is denned by the coefficient function, Agj 
l-C Bj (a s ). 



and observables. In particular, the four-loop APT expan- 
sion for the perturbative part Aej(Q 2 ) is given by 



(6) 



fc<4 



Here the coefficients Ck are the same as in Eq. (j3|), and the 
functions Ak(Q 2 ) are defined through the spectral func- 
tions Qk(c) = Im [a* (— a — ie)] by the spectral integral 



A(Q 2 ) = - [do-JL^ 

7T J O 





Q 2 



(7) 



Note, the first function, Ai(Q 2 ), is the analytic coupling, 
q:apt(Q 2 ) = A\(Q 2 ) . At large momentum transfers, all 
the functions Ak{Q 2 ) become proportional to the fc-th 
power of the usual perturbative coupling [a s (Q 2 )] k and 
the expansion © reduces to the power series ([3]). How- 
ever, at small enough Q < 1 — 2 GeV the properties of 
the non-power expansion ([6]) become considerably differ- 
ent from the PT power series ^ (see, e.g., Ref. Q for 
details) . 

2.2. The Q 2 -dependence 

Now we analyze the Q 2 -dependence of the BSR in the 
framework of both the PT and APT approaches in differ- 
ent orders (NLO, N 2 LO and N 3 LO) of the perturbative 
expansions ([3]) and ([6]), respectively. As a normalization 
point, we use the most accurate a s - value at Q = Mz, 
a s (M z ) = 0.1184 ± 0.0007 @, El]. In order to take into 
account flavor thresholds, we apply the matching condi- 
tions for the values of a s (Q 2 ) which are rather nontrivial 
in higher PT ordersjsee Refs. [H, [H ii|). Following 
to analysis in Ref. [23|, our matched calculation for the 
four-loop MS-coupling gives A ( ™' =3) = 336 ± 10 MeV. 
Note, we obtain practically the same results, but with 
larger errors, if we choose the pseudo-observable value 
R(Mz) = 1.03904± 0.00087 as a normalization point [H, 
which leads to the four-loop running coupling equal to 
a s {M z ) = 0.1190 ±0.0026. 

In Fig. [TJ we illustrate the behavior of the perturba- 
tive part of the BSR in different orders in a s in both PT 
and APT approaches. The APT curves in different orders 
(NLO, N 2 LO and N 3 LO) practically (at about 1 % accu- 
racy) coincide with each other, so we represent the APT 
result by a single dash-dotted line in Fig.Q] For complete- 
ness, we also show here the combined SLAC and JLab data 
on Y\~ n (Q 2 ) used in our analysis. The SLAC data points 
HI are denoted by squares, the JLab CLAS Hall A 2002 
data - by downward pointing triangles, the JLab CLAS 
Hall B 2003 data - by diamonds @, and the most recent 
JLab data [l[ - by circles. The horizontal dotted line rep- 
resents the limiting value T\~ n (Q 2 — > oo) = <?a/6. 

One can see that at Q 2 > 0.7 GeV 2 the four-loop ap- 
proximation describes the data quite well. Moreover, the 
corresponding curve passes close to the central values of 
several data points, although the experimental accuracy 
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Figure 1: Perturbative part of the BSR as a function of the momen- 
tum transfer squared Q 2 in different orders in both the APT and 
standard PT approaches against the combined set of the Jefferson 
Lab 0, and SLAC [H data. 

(which is of the same order as both the three- and four- 
loop contributions) does not allow one to make a definite 
choice between four- and three-loop approximations. 

At the same time, at Q 2 < 0.7 GeV 2 the four-loop 
approximation describes the data equally bad as the three- 
and two-loop ones. This a signal of the necessity to account 
for HT contributions, and it will be strongly dependent on 
the order of PT used for its extraction |4[ . 

This changes when APT is applied and the higher-loop 
stability is achieved. This is a well-known feature of APT 
free from unphysical singularities. At the same time, the 
deviation of APT curve from the data shows for necessity 
of the HT contribution which in this case is quite stable 

This situation may be considered as an indication of 
the transition of PT series to the asymptotic regime (while 
APT series remains convergent) for Q 2 ~ 0.7 GeV 2 . Let 
us explore this possibility in more detail. 

2.3. Convergence of the PT and APT expansions 

Clearly, at low Q 2 a value of the strong coupling is 
quite large, questioning the convergence of perturbative 
QCD series. The PT power series truncated after four- 
loop order (c.f. Eq. ([3])) reads 

AbTM = 0.3183 a s + 0.3631a 2 
+ 0.6520af+ 1.804 a* = ^ k{ot s ), (8) 

where Si is the i-th term. The qualitative resemblance of 
the coefficients pattern to the factorial growth did not es- 
cape our attention although the more definite statements, 
if possible, would require much more efforts. This observa- 
tion allows one to estimate the value of a s ~ 1 /3 providing 
a similar magnitude of three- and four- loop contributions 
to the BSR. 




Q 2 (GeV 2 ) 

Figure 2: The Q 2 -dependence of the relative contributions at the 
four-loop level in the PT approach. Four-loop PT order overshoots 
the three-loop one at Q 2 < 2 GeV 2 , so it does not improve the 
accuracy of the PT series compared to the three-loop one. 

To test that, we present in Fig.[2]the relative contribu- 
tions of separate terms in the four-loop expansion ([8]) 

N. i (g 2 ) = <5 l (Q 2 )/A Bj (g 2 ). (9) 
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Figure 3: The Q 2 -dependence of the relative contributions of the 
perturbative expansion terms in Eq. in the APT approach. Third 
and fourth order contributions amount to less then 5 % total, so 
the NLO APT approximation is sufficient for description of the low 
energy JLab data at the current level of experimental accuracy. 

As it is seen from Fig. [3J in the region Q 2 < 1 GeV 2 the 
dominant contribution to the pQCD correction Abj(Q 2 ) 
comes from the four- loop term ~ af. Moreover, its rel- 
ative contribution increases with decreasing Q 2 . In the 
region Q 2 > 2 GeV 2 the situation changes - the major 
contribution comes from one- and two-loop orders there. 
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Analogous curves for the APT series given by Eq. (j6|) are 
presented in Fig. [3] 

Figures [5] and [3] demonstrate the essential difference 
between the PT and APT cases, namely, the APT expan- 
sion obeys much better convergence than the PT one. In 
the APT case, the higher order contributions are stable at 
all Q 2 values, and the one-loop contribution gives about 
70 %, two- loop - 20 %, three- loop - not exceeds 5%, and 
four- loop - up to I %. 

One can see that the four-loop PT correction becomes 
equal to the three-loop one at Q 2 = 2 GeV 2 and noticeably 
overshooting it (note that the slopes of these contributions 
are quite close in the relatively wide Q 2 region) for Q 2 ~ 
1 GeV 2 which may be considered as an extra argument 
supporting an asymptotic character of the PT series in 
this region. 

In the APT case, the contribution of the higher loop 
corrections is not so large as in the PT one. The four- 
loop order in APT can be important, in principle, if the 
theoretical accuracy to better than 1 % will be required. 

2.4- The [i-scale dependence 

As it is known, any observable obtained to all orders in 
pQCD expansion should be independent of the renormal- 
isation scale /j,, but in any truncated-order perturbative 
series the cancelation is not perfect, such that the pQCD 
predictions depend on the choice of the /i-scale (for a re- 
view see, e.g., Ref. 

In order to estimate this dependence of rf~ n on the 
unphysical renormalization-scale parameter /i, we use the 
four- loop expression for the coefficient function Cej it^/Q 2 ) 
recently published in Ref. [6]. One commonly introduce 
the dimensionless parameter x^ (/i 2 = x^Q 2 ), which we 
have chosen to change within the interval x^ = 0.5 4- 2 
(see, for example, the analysis in Ref. [3]), and compare 
the /i-scale ambiguities between the three- and four-loop 
PT series. 

In Fig.21 the perturbative part of the BSR is plotted as 
a function of Q 2 in three- and four- loop orders of PT series 
corresponding to x^ in the interval 0.5 2. The width of 
the arising strip for the four-loop approximation is similar 
to the one for the three-loop approximation in the highest 
JLab region Q 2 ~ 3 GeV 2 Q so these approximations pro- 
vide the description of the data with comparable accuracy, 
as discussed above. Thus, the four-loop result does not 
improve the data description noticeably in the low-energy 
domain. 

At the same time, for Q 2 < 1 GeV 2 , where PT does 
not allow the description of the data, the inclusion of the 
four-loop contribution leads to a stronger /i-dependence. 

These observations provide yet other arguments sup- 
porting the mentioned transition to asymptotic PT series 
at Q 2 ~ 1 GeV 2 . 



2 One can find that an account for four- loop contribution leads to 
a decrease of the /^-dependence if Q 2 > 5 GeV 2 which is currently 
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Figure 4: The /i-scale ambiguities for the perturbative part of the 
BSR versus Q 2 for three- (shaded region between dash-dot-dotted 
and dashed lines) and four-loop (shaded region between short-dashed 
and dash-dotted lines) orders of PT series corresponding to changing 
of Xfj, in the interval 0.5 -f- 2. These two regions have similar widths 
and only slightly shifted w.r.t. each other, so the differences between 
three- and four-loop results are within the experimental error bars. 
Hence, in the common PT case, the N 3 LO approximation does not 
improve the data description compared to the N 2 LO one (see also 
Fig.0. 

3. Higher twists contribution 

3.1. The results of the fit 

Now, using expression @ fitted to the above men- 
tioned experimental data [2| we extract the coefficients 
/Z2i of the higher twist OPE corrections. The minimal bor- 
ders of fitting domains in Q 2 are settled from the ad hoc 
restriction x 2 < 1 and monotonous behavior of the result- 
ing fitted curves. 

Previously, a detailed higher-twist analysis of the two- 
and three-loop expansions in powers of a s was performed 
in Refs. Q. Now, we extend the analysis up to an order 
~ at- 

In Figs. [5] and [5] we present the results of I- and 3- 
parametric fits in various orders of PT and APT. The 
corresponding fit results for higher twist terms, extracted 
in different orders of PT and APT, are given in Table [T] 
(all numerical results are normalized to the corresponding 
powers of the nucleon mass M). From these figures and 
Table [T] one can see that APT allows one to move down to 
Q 2 ~ 0.1 GeV 2 in description of the experimental data 
At the same time, in the framework of the standard PT the 
lower border shifts up to higher Q 2 scales when increasing 
the order of PT expansion. This is caused by extra un- 
physical singularities in the higher-loop strong coupling. 



outside the JLab kinematical range, but will be accessible by JLab 
after the scheduled upgrade. 
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Table 1: Results of higher twist extraction from the JLab data on 




BSR in various (NLO, N^LO. 
APT. 



INTLO) orders of PT and all orders of 



Q 2 (GeV 2 ^ 



Method 


O 2 

^ mm i 






Ms /M° 


The best /M-fit results 


PT NLO 


0.5 


-0.028(5) 






PT N 2 LO 


0.66 


-0.014(7) 






PT N 3 LO 


0.71 


0.006(9) 






APT 


0.47 


-0.050(4) 






The best ^4,6,8-fit results 


PT NLO 


0.27 


-0.03(1) 


-0.01(1) 


0.008(4) 


PT N 2 LO 


0.34 


0.01(2) 


-0.06(4) 


0.04(2) 


PT N 3 LO 


0.47 


0.05(4) 


-0.2(1) 


0.12(6) 


APT 


0.08 


-0.061(4) 


0.009(1) 


-0.0004(1) 



Figure 5: The one-parametric /M-fits of the BSR JLab data in var- 
ious (NLO, N 2 LO, N 3 LO) orders of the PT and the all-order APT 
expansions. In the PT case, the four-loop result does not improve 
the data description compared to the three-loop one. In the APT 
case, the NLO approximation is sufficient due to higher-loop stability 
of the APT expansion (see also Fig- [3j . 
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Figure 6: The three-parametric /^4,6,8-fits of the BSR JLab data in 
various (NLO, N 2 LO, N 3 LO) orders of the PT and the all-order APT 
expansions. 

3.2. Sensitivity of the higher twists to Aqcd variations 

In the above analysis, we normalized a s at the Z-boson 
mass scale and then fixed the value of the A parameter 
separately in each order in a s approximation (it was suffi- 
cient for understanding the role of the fourth order in the 
PT/APT perturbative series). However, the correspond- 
ing values of the A parameter extracted in this way may 
be different from ones obtained in the direct QCD analysis 
of the experimental data on the moments of the structure 
functions (see, e.g., Ref. [26[). Having this in mind, we in- 
vestigate additionally the sensitivity of the extracted val- 
ues of the higher twist term ^,4 to the QCD scale parameter 
A in various orders of PT. In the framework of APT, the 
sensitivity of /14 to the A parameter is weak, and it does 
not depend on the order of the loop expansion. Corre- 



spondingly, the values of the higher twist coefficients turn 
out to be considerably more precise than those extracted 
in the PT approach (see also Table [T]). 




_q Qg I . [7^>0Oci\V:. , .llllll hi 

' 0.25 0.30 0.35 0.40 

A (GeV) 

Figure 7: Value of the higher twist coefficient 114 extracted from the 
JLab data using the PT at different orders at Q 2 mm = 0.66 GeV 2 
with error bands. Vertical lines denote the corresponding uncer- 
tainty ranges in A-parameter. The ranges corresponding to N 2 LO 
and N 3 LO approximations have similar sizes and overlap with each 
other, so the four-loop result does not improve the stability w.r.t. A 
variations compared to the three-loop one. 

In Fig. [7]we show values of the coefficient ^4 extracted 
from the JLab data using two-, three- and four-loop PT at 
Qmin = 0-66 GeV 2 vs the parameter A. One can see that 
the PT does not lead to a stable result for extracted \x\ 
value with respect to A variations. The extracted higher 
twist coefficient /Z4 changes quite strongly between differ- 
ent orders of the PT expansion. And it happens in both 
in absolute value and sign, namely, at A > 320 MeV the 
higher twist coefficient becomes positive in the four-loop 
PT order. This sensitivity of the higher twist term /14 to 
variations of the A becomes stronger at higher PT orders. 

On the other hand, these data tell us that the abso- 
lute value of fj,4 decreases with the order of PT and just 
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at four-loop order becomes compatible to zero. This may 
be considered as a manifestation of duality between higher 
orders of PT and HT (see Ref. Q and references therein). 
Moreover, when PT series manifests the asymptotic be- 
havior (i.e. becomes most close to exact result), the HT 
(which may be considered as a contribution completing the 
PT series) can be reduced to zero. 

4. Summary and Conclusion 

In this work, we performed the QCD analysis of the 
precise low energy JLab data on the BSR in the N 3 LO PT 
order and extracted the OPE higher twist terms using the 
four-loop expression for the QCD correction to the Bjorken 
integral Aej published recently in Ref. @. 

Our main observations are: 

i) The four-loop approximation provides good descrip- 
tion of the data for the highest JLab Q 2 ~ 3 GeV 2 . For 
several data points there is an impression that the four- 
loop approximation is better than the three-loop one. At 
the same time, the order of magnitude of both these con- 
tributions is the same as an experimental error, so a more 
precise statement can hardly be made. 

ii) For lower Q 2 < 0.7 GeV 2 the four-loop PT contribu- 
tion does not help to describe the data. Meanwhile, as it 
was shown earlier [J] , the APT application leads to higher 
loops stability of the HT extraction. In turn, this results in 
accurate data description down to Q 2 ~ 0.1 GeV 2 always 
at the two-loop APT level (see Fig.©. 

iii) The magnitude of HT decreases with an order of 
PT and becomes compatible to zero at the four-loop level. 

Our concluding impression is that all these features 
may indicate that the asymptotic nature of the QCD PT 
series is revealed at the four- loop level at Q 2 ~ lGeV 2 . 
This conjecture is confirmed by the analysis of relative 
contributions of various PT terms, as well as by that of 
unphysical //-dependence. 
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